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SYMBOLS ' 

Lower Case Latin 

a semi-diameter perpendicular to Z-axis in ellipsoid 
b semi-diameter parallel to Z-axis in ellipsoid 
f distributed loads in local coordinates 
h shell thickness; face sheet thickness 

i index: beginning edge of shell segment; independent joint of 

kinematic link; subscript "inside" ' - 

j index: ending edge of shell segment ; dependent joint of 

kinematic link _ . . ... 

n index on harmonic . . 

o subscript "outside" 

r radius 

s index of segment; coordinate in cylinder or cone 
t core thickness in sandwich shell 

V ‘ • ,r , 

w normal deflection, positive inward 

Upper Case Latin 


C stiffness eccentricity parameter's; offset distance in ogive, ellipse 

D bending stiffness parameters 

2 ’ 

E Young’s modulus (lb/in ) 

F lineal force (lb/in) 

G shear modulus (lb/in ) 

K extensional stiffness parameters 

M ; bending moment on shell- (in-lb/in) 

N membrane stress resultant (lb/in) 

N assumed membrane stress resultant (lb/in) 

Q transverse shear stress resultant (lb/in) 

R radius; "global" coordinate, positive radially outward. 
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SYMBOLS (continued) 


T temperature; "global" coordinate , tangential 
X Cartesian coordinate, 9 = 0 at X-axis . 

Y Cartesian coordinate 

Z Cartesian and "global" coordinate, coincides with axis of revolution 

Greek. 


a angle between rotated coordinates 

P ratio of semi-diameter parallel to Z-axis in ellipsoid to semi- 
diameter perpendicular to Z-axis 

7 shear strain; non-linear parameter; angle of inclination of 
kinematic link 

C normal coordinate, positive inward r - 

9 circumferential angular coordinate (rad) • • 

A shell parameter 

v Poisson's ratio 

2 

a normal stress (lb/in ) 

2 

x shear stress (lb/in ) 

<f> meridional angular coordinate (rad) 

to rotational displacement (rad) 

'rotational displacement in "global" coordinates (rad) 

f 

A displacements in fixed or "global" coordinates 
A segment length parameter 

Miscellaneous 

eq equivalent 

s4> sin p 

cp cos <p 

Other symbols are defined in the text where used. 
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SECTION 1 


PROGRAM CAPABILITY 


The use of an accurate shell theory to analyze structural shell problems 
usually involves complex mathematics and numerical techniques, which are 
nearly impossible to treat without the aid of automated procedures. On this 
basis, digital computer programs based upon the Love-Reissner first order 
shell theory have 'been developed. These programs can perform a stability or 
vibration analysis of ortho tropic thin shells of revolution, subjected to 
axisymmetric distributed loading or concentrated line loads, as well as 
thermal strains (Reference l). Furthermore, a shell with arbitrary boundary 
conditions, underloads which vary arbitrarily with position and under a 
temperature variation through the thickness, is tractable with these .programs. 
The shell can consist of any combination of the following geometric shapes: 

1) Ellipsoidal - spherical (offset from the axis of revolution allowed) 

2) Ogival - toroidal 

3 ) Modified ellipse shape 

4) Conical - circular plate 

5 ) Cylindrical 

6 ) General point input geometry 

7 ) Dummy geometry slot to be filled by the user 

8 ) Discrete ring 

The shell wall crossection can be a sheet, sandwich, or reinforced sheet or 
sandwich. The reinforcement can consist of rings and/or stringers, a waffle 
construction rotated at any angle to the principal coordinates, or an isogrid 
construction. General stiffness input options are also available. The rein- 
forcement material properties can differ from those of the main shell, and a 
temperature variation can cause different properties in the two face sheets of 
a sandwich shell. 
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The basic approach to the problem (Reference l) is to cut the structure into 
several shell regions. These regions need to be singly-connected shells, 
and can only have line loads applied at their end points.' There are no 
restrictions on geometry, or uniform or thermal loads. The regions are 
further subdivided into several shell segments, each being free to have its 
own geometric shape, provided that the shape falls into one of the categories 
mentioned above. 

Stiffness matrices obtained for each segment, are coupled by standard matrix 
methods to obtain region stiffnesses, which, after being reduced in size, are 
in turn coupled to form the total shell structure under analysis. Currently, 
the UNI VAC 1108 computer programs are sized to handle a structure composed of 
up to 29 segments in each of 29 regions arbitrarily connected to each other. 
There is a limitation on the size of a shell segment, which is a consequence 
of the demand that boundary disturbances be felt throughout the segment. This 
limitation is mathematically described in Section 2 (pages 2-31 to 2-33) as a 
length parameter. This parameter, however, is not reliable near the apex of 
any shell shape (<J> = 0), and the segments needed in this region are actually 
much smaller than predicted by the parameter. A mathematical singularity 
occurs at the apex where r Q (the radius of revolution) becomes zero. It is 
this singularity which prevents the length parameter from being meaningful 
near the apex. Furthermore, the point (<j> = 0) is not ah acceptable input 
point of the program (except for the torus-ogive and offset ellipsoid), 
although any point outside a circle of infinitesimal radius is satisfactory. 

There is a considerable latitude in what can be done within each shell seg- 
ment. The thickness of any segment can be symmetrically tapered and it can 
contain up to 14- points of discontinuity, provided that the segment center- 
line remains continuous and describable by a single shell geometry. A temp- 
erature distribution through the thickness can be specified at three points 
in a homogeneous shell, and 4- points in a shell of rigid core sandwich con- 
struction. The distribution is considered to be linear between these points. 
Thus, it is possible to approximate temperature distributions other than 
linear distributions. In the event of physically discontinuous shell center- 
lines, a kinematic link is available for use in the analysis. The link 
relates displacements across the discontinuity. This link may be used between 
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regions, and between segments within a region. Discrete offset rings are 
also available for use within or between regions. 
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SECTION 2 


INPUT INFORMATION 


2.1 GENERAL NOTES 


The preceding section provides some insight into the capability of the pro- 
grams, and the potential that they might‘ v have for future use. If the pro- 
grams are applied judiciously they can be extremely powerful tools. The 
mechanics of applying them should be clearly understood. With this in mind, 
the remaining section should be studied carefully. 

The required input data may be subdivided into three main parts, namely: 
geometric, topological (or coupling orientation) and joint data (degree of 
freedom description for each joint component). Each segment requires its 
own geometric configuration and numerical integration control. 

The output consists of stiffness coefficients for each shell segment and the 
actual symmetry of the coefficients is presented in a convenient form for 
a check on the accuracy of the integration through the segment. Region 
stiffnesses and their symmetry checks are also provided. Eigenvalues and 
final eigenvectors are printed out for each shell segment at intervals 
along the segment as specified by the user of the program. The output will 
be further discussed in Section. 3. 

The present program size is described in the table below. 

UNIVAC 1108 Table of Program Sizing 


I. 

Segments- per region: 

29 

II. 

Segment joints per region: 

30 

III. 

Regions: 

29 

IV. 

Region joints: 

30 
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UNIVAC 1108 Table of Program Sizing (continued) 


V. Number of points available 
per segment for specifying 
geometric or load data: 

VI. Number of points available - 
through the thickness for 
specifying temperature data: 

VII. Geometries: 


VIII. Wall cross-section options: 


IX. > Number of material property 
tables per submission: 


X- 


XI. 


-Number of points per 
material, property table : 


A combined fixed and 
variable loading option 
is available. 

XII. Orthotropy options: 


30 


ellipsoid, sphere, offset ellipsoid, 
modified ellipsoid, ogive, toroid, 
cone, annular plate, cylinder, general 
geometry, ring, elastic support, dummy 
geometry. ' r - ; - 

single sheet, equal face sheet sand- 
wich, unequal face sheet sandwich, 
eccentric reinforcement (rings, 
stringers or. both:) ,w.af f lerxeinf orce rv' 
ment rotated at an arbitrary angle to 
coordinate axes, isogrid reinforcement, 
arbitrary stiffness 1 input: ■ l ’ K •• r ' ? - 


10 


10 


isotropic or orthotropic sheet, iso^ - 
tropic or orthotropic sandwich, iso- 
tropic or orthotropic sandwich with 
different face sheet .properties caused 
by thermal gradients, isotropic or 
orthotropic sheet or sandwich reinforced 
by different property rings or different 
property stringers or both, isotropic or 
orthotropic sheet or sandwich reinforced 
by a different property waffle system 
rotated by an angle (3 to coordinate axes, 
isotropic or orthotropic sheet or sand- 
wich reinforced by a different property 
general isogrid reinforcing system, 
arbitrary stiffness input options used to 
describe , other configurations. 
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UNIVAC 1108 Table of Program Sizing (continued) 


XIII. 

Number of rings per region: 

28 


XIV. 

Total number of rings at 
region joints in structure: 

28 



Figure 2-1 shows the detailed option flow chart for the present programs. 
GENERAL MOTES - Idealizations ( . 

Before discussing the. specific card input order, it would be advantageous 
to introduce some general guidelines in the area of idealizations and 
topology. In many computer programs there is such an abundance of numerical 
computation, that minimizing numerical roundoff errors becomes as important 
as getting the final answers . In some cases the engineer can aid the 
program in this effort throught the use of judicious idealizations. Such 
a possibility exists in the STARS-2 programs, since many internal operations 
are 'involved with building and inverting stiffness matrices. The object 
of the user therefore, should be to help the computer by avoiding the 
creation of ill-conditioned matrices at any step (see Reference 2). 
Physically, the way to achieve this end is to have all the segment stiff- 
ness matrices of the same order of magnitude. This will in turn produce 
region stiffness matrices which are of similar orders of magnitude, and 
minimize possible ill-conditioning in the total structure matrices. The 
.user can help to achieve this end by sizing his segments in such a way 
so that no short stiff segment is contained alone in a region with all 
other long flexible segments, or that no region comprised of all short 
stiff segments exists in a structure whose other regions contain only 
long flexible segments. No accurate measure' can be given on the relative 

stiffness or flexibility of segments allowed, and thus the best check is 
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Figure 2-1 Program. Option Flow Chart 


2-4, 2-5 


































to see if a structure is in equilibrioum under the applied loading. It 
also must be kept in mind, that if an idealization has provided useful 
results for axisymmetric (zeroth harmonic) loading, it need not neces- 
sarily be a good idealization for other harmonics, but if good results 
are obtained for the first harmonic, then other harmonics should also 
encounter no adverse behavior. The symmetry checks of segment and 
region stiffness matrices are useful for many reasons, but will not 
necessarily alert a user to ill-conditioning. 

In the use of regions, one other type of accident must be avoided. This 
is the creation of a single region structure with both ends fixed, where- 
in no suitable boundary condition matrix can be formed. Thus, in the 
use of region idealizations, which are less physically meaningful to a 
user than pure segment idealizations, care should be taken so that all 
boundary conditions are not zeroed out. To avoid this problem, and to 
minimize program running time, it is best to maximize the number of 
regions in a structure, and minimize the number of segments per region. 
Thus, in small problems, for best numerical efficiency, there should 
only be one segment per region. 

In the solution of the static axisymmetric problems, and the n = 0 
eigenvalue search, all torsional degrees of freedom should be removed. 

A second boundary condition matrix should be input for searching for 
eigenvalues in other (n f 0 ) harmonics. 
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GENERAL NOTES - Data debugging 

The STARS programs have been provided with special separate data debugging 
packages called SATELLITE programs. In order to be able to debug as much' 
of a given data deck as possible in one computer submission, the data is 
grouped by inserting special cards, termed "dash-separator cards", appro- 
priately. In order so that additional errors are not made by requiring 
insertion and removal of these cards, the STARS program has been coded to 
accept these dash-separator cards in the input. A dash-separator card is 
shown below: 



As can be seen, a minus symbol is inserted straight across the computer, card, 
from column 1 through column 80. 

Since the dash- separator cards subdivide the data deck, there exists the pos- 
sibility that a separated data block may be completely omitted (for example 
no kinematic links in a structure). - In this case one dash-separator card is 
also omitted . Under no circumstances can there exist two adjacent dash- 
separator. cards in a data deck. The SATELLITE programs are described in 
Reference 3* 
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DETAILED ORDER OF INPUT (See Figures 2-1, 2-2) 


GENERAL 

INTRODUCTORY cards 

Column 

Forma -1 

1. 

Title Card 




A. 

Alphameric title (submission description) 

i-64,: 

16A4 

2. 

Program Control Card 




A. 

Number of regions to be coupled 
(Max. = 29) 

1-2 

12 


B. 

Total number of segments 
(Max. = 29 x 29 = 84l) 

3-5 ■ 

13 


C. 

Number of Material Property Tables 
(Max. = 10) 

• 6-7 

12 


D. 

First harmonic in the range to be searched 
for stability or vibration eigenvalues. 

■ 8 “9 

12 


E. 

Last harmonic in the range to be searched 
for stability or vibration eigenvalues. 

10-11 

12 . 


F. 

Search harmonic increment . 

12-13 

12 

V 

G- 

Number of loading conditions for this problem 
(See Dynamic Option Card for Vibration 
problems.) 

14-15 

12 


In a stability analysis two types of loadings 
may be applied simultaneously., The first 
loading will be considered as being of 
constant fixed magnitude, while the second 
will be considered as being the variable 
(eigenvalue) loading. In this case the 
number of loading conditions is set to 2. 

If there is no fixed load applied, then 
the number of loading conditions is 1. 

H. Number of boundary conditions for this - 16-17 12 

problem. 

In some cases (such as shells containing an. 
apex) boundary conditions must be changed " ' 

for the first few harmonics. An option thus " 
exists of using up to 3 different boundary >• • 
condition . sets in an eigenvalue analysis. 

For example, if the first eigenvalue sought''' 
is for the n =. 6 harmonic, torsion degrees " 
of freedom should be removed. For eigenvalue 
searches in other harmonics torsional freedom 
should be retained by using a second set of 
boundary conditions. 
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. Column •. 

I. Prestress input clue 18-19 

If this clue is set to zero, the program 
will expect load' input, and will' calculate 
the prestress state by itself. In 
considering the prestress state it will 
include the effects of predeformation. 

If this clue is set to unity, the program 
will^ expect the prestress state to be input, . 
and will assume there is no predeformation. 

If there are two loading conditions (item 
G above), the program will expect two sets t ■' 

of input prestress states. ! 

The prestress' input clue set to unity may 
also be used to obtain estimates of buckling • 
under unsymmetric loading by inputing the 
"worst meridian" . stress distribution as an 
axisymmetric condition. This technique has 
been justified experimentally for cylinders 
under some loading conditions DO- 

J. Intermediate print' clue 20-21 

If this clue is set to unity, base stiffness 
and prestress stiffness matrices for segments, 
regions, and' structure will be printed for all 
estimates in all harmonics. To delete this 
intermediate print the clue should be set to 
zero. " . ” . 

K. Eigenvalue calculation clue . . 22-23 

If this clue is set to zero an eigenvalue matrix 
procedure will be used for buckling load or ' 
vibration frequency calculations Cl]. If this 
clue is set to unity a determinant evaluation 
procedure will be used Cl]. Since the 
determinant evaluation procedure is much 
longer and can skip roots based on user input, 
it is not recommended for practical analysis. 

L. Convergence Tolerance , 24-30 

If- this is .left blank .and an eigenvalue procedure 
is used for calculations, the iterations will' be' 
carried on until successive estimates of the 
eigenvalue differ by less than 1.0$. Otherwise 
the user can input his own convergence tolerence . 

(in decimal form, not $).' 


Format 

' 12 


12 


12 


F7.0 
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Figure 2-2 Data Sequence 


2-16 







Column 


Format 


M. Upper cutoff value 31-37 

This input is mainly relevant for stability 
calculations although it will be used for 
vibrations if input. If the user wishes to 
speed his computer search time by disregarding 
any eigenvalues higher than a specific ratio 
on his input prestress state, he may input 
this ratio here. Thus if 2.5 is input, all '• 

•“ harmonic calculations where the minimum 
.. eigenvalue is found to be greater than 2.5x 
• (input prestress state) will be aborted. The 
~ cutoff value will not be u£ed if this field 
,_is left blank. 

N. Determinant evaluation largest increment 38-44 

' Input here and item 0 will be disregarded if 
-the eigenvalue calculation option is used 
(item K above). 

This is the initial load step (input as a 
decimal, used as a percent) used for the l 
loading variation in the determinant 
calculation option. If this field is left 
blank the program will use 5 $. 

O. Determinant evaluation smallest increment 45-51 

multiplier 

This is the multiplier for the initial load step 
which will determine the smallest load step. 

The multiplier can start at 1.0 and should 
decrease in a factor of 10 to (for example) .1, 
or .01, .001 etc. Thus for an input of .05 
above ( item N ) and . 01- here , the load will 
first be changed by 5#-, then .5$,. and finally 
.05$ (=5$x.0l) as convergence increases. If 
this field is left blank the program will use 
1.0 making the smallest step = to the initial step. 

P. Harmonic Mass Multiplier (for vibrations only ) 52-58 

To use the harmonic mass multiplier see 
Reference 1 and the Dynamic Option Card. This 
entry is not -necessary, for a free or prestressed 
vibration problem, or a free or prestressed 
critical speed problem. 

Q. Frequency or Angular Velocity in rad. /sec. (for 59-72 
vibrations only ) 

The assumed value of the frequency of 
interest (see Dynamic Option Card). 


F7.0 


F7.0 


F7.0 


F7.0 


E14.7 
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Column 


Format 


3. Dynamic Option Card (omitted for stability ) 

A. Dynamic Analysis Description Clue . 1-4 A4 , 

There are six possible dynamic.- analysis options, ' 1 -• i 

each defined by a code word. If the run is to be ' ‘ 

a free vibration analysis (unloaded shell), the • - • :-.c ■' 

code word is FREV and item G on card 2 should be . • i ’ - 

zero. 

If the run is to be a vibration analysis of a pre- . ; . 

stressed shell, the code word is to be VPRE and 
item G on card 2 should be unity. If the run is to 
be a critical rotation speed analysis, the code • i " 
word is CRSP and item G on card 2 should be unity... 

If the run is to be a critical rotation speed 
analysis of a prestressed shell, the code word 
is PCRS and item G on card 2 should be two. 

If the’ run is to be a critical rotation speed 
analysis including mass perturbation (see Ref . 1), 
the code word is i CRSR and item G on card 2 should . 
be unity. ! 

If the run is to be a critical" rotation speed ‘ 
analysis of a prestressed shell including mass 
perturbation (see Ref. 1), the code word is PCSR 
and item G on card 2 should be two. 

B. Eigenvector Clue 5-10 16 

If this clue is zero the lowest two eigenvectors 
will be sought and printed for each harmonic. If _ . 

this clue is unity the two eigenvectors corres- 
ponding to the frequencies closest to the estim- • .. 
ated frequency will be sought and printed for- 
each harmonic. 

C. Rigid Body Mode Clue 11-15 15 

The n - 0, 1 harmonics contain rigid body modes 
for a totally unsupported (free) shell. If an 
analysis for other modes in these harmonics is 
required for this , type of shell, the clue is -set 
to unity. In this case the total harmonic search 1 
can only encompass n = 0,1. The first. two eigen- . 
values in the output should then approximate zero, 
and the third will be the first to correspond:!, to a 
nonrigid body mode. In all other cases (including 
a free shell investigated for harmonics n>2), 
the clue is zero. 

MATERIAL PROPERTY TABLES (Max. = 10 sets ) 

As many sets of these cards are used ( 10) as there are different 

material property segments in the structure , to be analyzed. These . 
tables will be used to obtain the thermal variation of material 
properties if thermal loadings exist. Thus the range of temperature 
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- Column 

in this table should be greater than that of the thermal , 

loads. If no thermal loads exist, the values given in the 

first column of this table will be used, and the rest of the 

table can be left blank. . If. there are thermal, loads,, the 

range of the table is to be considered as that between the , „ 

second and tenth columns. • • 

1. Identification Card • ‘ • ■ < . 

A. Material Title (Alphameric) , ' - 1-^ 

Any name can be made up as long as it. is con- ■- i . . . ; . 

sistently used on the segment cards to which • 
it refers. The same name cannot appear on 
more than one (l) table. 

. i \ i 

B. Type of Table ....... • , . . , 11-14. . 

One of several possible alphameric’ clues is 

written here. ’These clues serve to' size the •’ 

number of cards in the property table, and 

define which properties belong on which card. . 

The possible clues are: 

ISOT ' ’ • - 1 • 

ORTH 

STIF 

Their definitions are provided in Item 2 below. ... 

2. Material Property Cards 

The material property cards below are given depending 
upon which table type clue is used* If the table 
type clue is "ISOT" (isotropic table): 

A. Temperature values (5 values per card; 2 cards) 

These are the temperatures at which the values 
of material properties will be given. The first 
value in the table must always be the room or' 
stress-free temperature, ' since ’the material- 
properties in- only- the first column of the table 
will be used- in an analysis involving no thermal 
load. The rvalues of temperature in table col- 
umns 2 through; 10 must be in ' algebraically 
increasing order. 

B. Values of Young’s Modulus at the given temperatures. 

(5 values per card; 2 cards) 

C. Values of Poisson's Ratio at the given temperatures. 

(5 values per card; 2 cards) 

D. Values of the thermal coefficient of expansion 
at tlie given temperatures,. 

(5 valups per card; 2 cards) 


Format 


Ak 

Ak 


5E14 . 7 


5El4 . 7 
’ ' 5E14 . 7 

5E14.7 


2-19 



MATERIAL PROPERTY TABLES (continued) Column Format 


If the table type clue is "ORTH" (orthotropic table): 


A. 

Temperature values - 
(5 values per car’d; 2 cards) 

These are the temperatures_.at' / which, the values of 
material properties will be given. 

5E14. 

B. 

Values of Young's Modulus in the 9 direction (Eg) 
at the given temperatures. - 

(5 values per card; 2 cards) 

5E14 . 

C. 

Values of Young's Modulus in the <{> direction (E ). • 
at the given temperatures. 

(5 values per card; 2 cards) 

5El4. 

D. 

Values. of the Poisson' s, Ratio v at the given 
temperatures. ' • 

(5 values per card; 2 cards) 

5El4 . 

E. 

Values of the thermal coefficient of expansion in 
the 9 direction (a Q ) at the given temperatures. 

(5 values per ( card; 2 cards) , 

5El4 . 

F. 

Values of the thermal coefficient of expansion in 

the <f> direction (a x ) at the given temperatures. 

9 

(5 values per card; 2 cards) 

5El4. 

G. Values of the Shear Modulus G, „ at the given 
temperatures. * 

(5 values per card; 2 cards) 

If the table type is "STIF" (table to be used for 
reinforced shells): 

5E14 

A.-G. 

The values in these locations are the same as 
those above for the "ORTH" clue case, and refer 
to the basic shell. 

5El4 

H. 

Values of ring Young's Modulus (E^) at the ■ 
given temperatures. 

(5 values per card; 2 cards) . 

5Ei4 

I. 

Values of stringer Young's Modulus (E ) at the' 
given temperatures. 1 

(5 values per card; 2 cards) ! 

5El4 

J. 

Values of ring thermal coefficient of expansion 
(a^) at the given temperatures. 

(5 values per card; 2 cards) 

5Ei4 

K. 

Values of stringer thermal coefficient of 

expansion (a ) at the given temperatures, 
s 

(5 values per card; 2 cards) 

5El4 
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Column 


Format 


. Note : In a rotated waffle or isogrid construction, 

items H and I, and J and K, refer to the. grid 
directions and are respectively identical. 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus in 1-80 

(See General; Notes - Data Debugging) 



REGION INTRODUCTORY CARPS ' . .. 

These two cards are placed at the beginning of each region 
data information. Each region contains the following data : 
set (see Figure 2-2): a) Two region introductory cards; 

b) data cards for each segment within the region; 

c) ring cards describing the discrete rings within the 
region, if any; and d) kinematic link cards' describing 
the kinematic links within the region, if any. 

1. Identification Card 

A. Number of segments within the region- (< 29 ) 

B. Number of kinematic links ‘ between segments 
within the region. 

C. Number of discrete rings between segments - 
within the region. 

D. Any alphameric information (region description) 

2. Topology Card (Coupling Orientation) 

A. Region Number 1-5 15 

Number of the region under consideration. 

B. Joint (i) ^ 6-10 15 

Joint associated with i (beginning) end of 

the region. 

C. Joint (j) 11-15 15 

Joint associated- with j . - (ending) end of 

the region. 

There is no coordinate flow in regions (unless 
1 region = 1 segment), such as that shown for 
segments in Figures 2-3 to 2-9* However, the 
start joint of a region must match with 1 in 
segment numbering, and the end joint must 
match with' the highest' segment joint number 
in the region (see Figure 2-12 and page 2- 5 I 4 .) . 

If 1 region = 1 segment the segment topology 

card will be a dummy card [ 112 ] (see page 

2-54) - ' / '. ’ ’ 


1-2 

12 

3-4 

12 

5-6 

12 

7 -TO 

i6a4 
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Figure 2-3- Typical Shell Segment 
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Figure 2-5. Modified Ellipsoid 
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Specify’ r^, C ;> 0 



Figure 2-6. Ogive 



z 



Figure 2-7. General Geometry 
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Figure 
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out 



Specify r Q 



-- 9 * Cylinder 
.2-29 



SEGMENT CARDS Column 

This sequence of cards is repeated- for each segment 
within the region. 

1. Identification Card 

A. Segment identification code (see Figures 2-4 to 2-9) 


1. Ellipsoidal or spherical shell Code = 11 1-2 

2. Modified ellipse shape Code = 12 1-2 

3. Ogival - Toroidal 'Code = 13 1-2 

4. General Geometry ' Code = l4A or 1-3 

(see Figure 2-7) 14 b 

5. Dummy geometry slot Code = 15 1-2 

( 4> coordinate) 

6. Conical - Circular Plate Code = 21 1-2 y 

. The plate is treated as a cone with zero angle.- - - 

7. .Cylindrical shell Code = 31 1-2 

B. Any alphameric information (segment description) 4-67 

2. "MAGIC" Control Card 

A. Interval at which final answers are to be printed l-l4 

( in radians or inches ) . 


The <j>-coordinate is defined for all geometric 

shapes except the cylinder, cone and plate, 

for which the s coordinate is used. Figures 

2-4 through 2-9 describe these coordinates - • ' - 

for each shape. 

B. Difference 15-28 

The value recommended depends upon the 
computer used. For eight figure accuracy 
computers it is 1.0 E-6; for the IBM 3&0 
it is 1.0 E-4. 

,C. Integration interval 29-42 . 

The Runge Kutta numerical integration procedure 
is substantially more accurate than finite dif- 
ferencing. An interval of (.,01 to. 03) x segment 
size (in radians or inches) is recommended in eigen- 
value analysis so as to be able to represent eigen- 
vectors with high wrinkling in the segment interior. 

(in using a 30 point segment table (see p.2-37)there 
should be at least 30 integration steps.) 

D. Delta . . 71-72 

For a fixed -step integration. Delta = 0. 

This card controls the Runge -Kutta numerical 
integration scheme. The suggested values above 
yield accurate results for a fixed-step integration 
method . 


j Format 


F2.0 

'F2. r 0 IJi - 

F2.0 

F2)0Al 

F2 . 0 " 

F2.0 , 

F2.0 

i6a4 

E14.1 


E14.1 


E14.1 


, F2.0 
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. Calculation .of Segment Length 


There is a restriction on the length of the shell segments. Physically, 
the restriction demands that boundary disturbances at one edge b’e 
distinctly felt at the other edge. This is a consequence of using a 

numerical integration procedure. Since the segment stiffness matrices 

C . ' ... 

must be symmetric, the calculations involved in obtaining each matrix element 

r\ " ’ . 

■must be such that a computer rpund off error never becomes prominent. Limit- 
ing the segment length insures satisfaction of this criterion.' This length 
is a function of both geometric shape and segment location -within a specific 
geometry. One of the limiting factors is that the ratio of the radii of 
revolution at the initial and final points' of a segment be greater than one 

hundredth and less than one hundred. Thus I <100 where: 

: — 100 Wj / 



This requires smaller segments than will normally be predicted by formula in 
the area of an apex. In addition, note that (<j> = 0) is not an acceptable 
input point (except for the torus-ogive or offset ellipsoid). 

For a cylinder, the segment length parameter, ■ 

A = (1 + y) 2 0 As 


should be held to about ^.0. In this expression, "7" is a non-linear param- 

. \ , . . . * - ; 
eter. For homogeneous shells:- 




It is zero for a linear problem.' 


"0" i‘s a measure of the rate of decay of a disturbance in the shell. 

: . > 

"As" is the meridional length, 
r - r_ sin <p 

O 2 ' 
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The values of 3 and As for various shell geometries are given below: 




For w • 0.3, A s k 



i 

Hcaofeneoua Cylinder 
1 

H 

T 

S-X^rt 1 

3.11<r.B) 2 . 

A. S 1 * 

, ; r) 

i 

Sandwich Cylinder - 

-Equal Face Sheets 

, k . Ul - 

■ ' f 0 2 (« ♦ 6bt ♦ 3t 2 ) 

A. i 3.11 £ r 0 2 (kh 2 . 6ht ♦ St 2 ! j. 

-L. 


. J .£ 

I* 

h 

T 


. * 

Sandwich Cylinder - 
Unequal Face 8heete 

- ju-. 2 ) r h i * 

r | £ 

(bj ♦ h o )'* • 12b 1 b () t(b 1 ♦ h 8 ♦ t) L r 0 J 

A. l 3.11 . ,J [<», * b 0 > ♦ lOhjb^lbj . b o . t>] 

* 



T ° 


. 


Approximate formulas can be obtained for near cylindrical regions of generally 
curved surfaces. The length parameter. 


A = (1 + y)* 5 A A<|> 

should be held to about 4.0. In this expression "y" has the same definition 
as in the cylinder case.- ' 

"A" is a measure of the irate of decay of a disturbance in the shell. 

As 1 

"A<f>" = — — is the angle intercepted by a meridional arc length "As". 
r l 

4 - 

The values of A and As for various shell geometries are given below: 




k 

. 2. F 1 






Sandwich Construction - 
Equal Face Sheets 


3(1 - v 2 )^ 1 * 



r 2 2 (Ab 2 ♦ «bt ♦ 

3t*) 

Sandwich Construction - 
Unequal Face Sheets 

l k - 

3(1 - 

r b i- b =r 

(b l 4 ho 11 ' 4 l2h i l, o tlh i ♦ h 0 ♦ t > L r 2 J 


For V • 0.3, A i b: 


3.ll(r,B) 
A. a &r 

(l * yP 


A 

2 


At 1 3.11 [ r 2 J (Ab 2 ♦ 6bt * 3t 2 )J 


t 


[b^-t] [ (b i ♦ b o )k ‘ “wk ♦ b „ ♦ *>] 
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I 


' _3 

)The minimum allowable segment length is 1 x 10 
j(inches or radians according to segment sizing). 

| ~ • Column Format 

h. Geometric Description Card 

A. Ellipsoid and sphere (Figure 2-4 ) 


4- 

1. Semi-axis perpendicular to Z-direction (a) 

2. Ratio of semi-axis in the Z-direction (b) 

i-ih 

15-28 

E14.1 

ei4.i 

to (a), P = (|) 

3. C = offset distance (±) (C = 0.0 if no offset) 

29-42 

E14.1 

B. 

Modified ellipse shape (Figure 2-5) 




1. Axis ratio coefficient (n) 

2. Semi-axis ‘perpendicular to Z-direction (a) 

1-14 

15-28 

' E14.1 
E14.1 

c.- 

'Ogive (Figure 2-6) 
1. r = radius 

1-14 

E14.1 


2. C = offset distance (±) 

15-28 ' 

E14.1 

D. 

General Geometry (Figure 2-7) 

1. Number of pairs of Z versus r Q points (<l4) 

1-2 

12 


2. Z versus r Q points, in pairs, starting 

3-72 

F10.0 


with Z value (7 values per card, including 
first card, for up to 4 cards total). 

1-70 or 

F10.0 


In the input table the first Z value is 

taken as Z. = 0, and furthermore Z. > 0 
.1 l 

(i = 2 - l4) (see Figure 2-7). 

E. Cone (Figure 2-8) 

1. Angle <f> in radians (for flat plate, <f> = 0). 1-1 4 E14.1 

Keep in mind that this Ms a constant for 
a given cone and should not be confused 
with the <f> on the MAGIC Control Card. 

{ ’ 

F. Cylinder (Figure 2-9) 

' ” 1. Radius ' " l-l4 “ El4.1 

4. Master Clue Card 

: This card contains a series of clues which de- 

| termine the program and table locations to be 
« used for the segment being described. For a 
master flow chart of clues and options in the 
program see Figure 2-1. 
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Column 

Note : In a two loading condition stability analysis 

only one condition can involve temperature. Thus the - , - 

thermal clue on this card refers to that load condition.. 

One load clue card in this case must contain a zero in " . 

column 1 (see capd set 8). 

A. Material Table Type Clue , 1-4 

This clue defines the type of material property 
table to be expected for the segment; This, as ' 
well as the following clue determines the mate'- 1 ■ 
rial properties that will be used in the struc- 
tural analysis for the segment. Thus -these two 
clues should match the two clues used' on the '•> ; ' 

identification card of the corresponding mate- • ' 

rial property table-. As mentioned before -on ;j - .- 

page2-19, the three possibilities are:' '■ 

; isot * . . ■ , . ’ 

ORTH " . '. . -. 

STTF . , _ ... • ' , 

B. Material Title ' ' ' ll-l4 - 

This name should be the same, as the name which, .. 

appears on the material property, table which 

contains the' properties to be utilized for 

this segment. ' . ’ .. ' , 

C. Sheet Clue * " •••••'•' -21-24' 

- i . • 

. ■ I . . I 

This clue informs the program as to what kind . 

of shell wall crossection to expect. If. the . , 

shell is of single sheet construction, the clue 
to be used is: SING. If the shell wall is an 

equal- size x f ace sheet sandwich, the clue to be 
used is: BfepUA . If the shell wall is a sand- 

wich but the face sheets are not equal) the 
clue to be used is: UN E Q,. < Finally, if- the - • 

shell segment is reinforced by -rings, stringers, a 
waffle, or an isogrid, the ■ clue to be used ' '■ 

is: BLAI-J . . - - - 

D. Reinforcement Clue 31-3^- 


Format 


A4 


A4 


A4 


a4 


This clue describes the type of reinforcement 
that is present on the .shell. If the shell is 
purely of single sheet or equal or unequal-size 
face sheet honeycomb construction (no reinforcing),., 
the clue to be used is: THIC. . If the reinforce- 

ment consists of rings or stringers or both, 
located along the coordinate axes (0 .and <p or s), 
three clues are possible .depending upon the 
basic shell wall construction. If the basic 
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Column 


Format 


wall construction is a single sheet, the clue to 
be used is ST11 . If the basic shell wall is an 
equal-size face sheet sandwich, the clue to be 
used is ST12. If the basic shell wall is a sand- 
wich but the face sheets are unequal, the clue to 
be used is ST13 » If the reinforcement consists 
of a waffle which is rotated at an arbitrary angle 
3 from the meridional axis the following three 
clues are possible depending upon the basic shell 
wall construction. If the basic wall construction- 
is a single sheet, the clue to be used is RWA1. 

If the basic shell wall is an equal-size face 
sheet sandwich, the clue to be used is RWA2. If 
the basic shell wall is a sandwich but the face 
sheets are unequal, the clue to be used is RWA3 . - 
If the reinforcement consists of an isogrid con- 
struction of general angle (3 from the meridional 
axis (normally P = 30°) the following three clues 
are possible depending upon the basic shell wall 
construction. If the basic wall construction is 
a single sheet, the clue to he used is ISG1 . If 
the basic shell wall is an equal-size face sheet 
sandwich, the clue to be used is ISG2. If the 
basic shell wall is a sandwich but the face sheets 
are unequal, the clue to be used is ISG3 *' Two 
other clues are available, namely ST10 and RWAF, 
if the user wishes to input his own stiffness 
constants. These constants may represent any 
wall construction as long as the basic Hooke's 
Laws used with the clues are appropriate to 
describe the construction to be considered. The 
Hooke's Laws used with these clues are given 
under the description of segment card set 6. 

Note : The reinforcement described in the seg- , 

ment cards' is closely spaced reinforcement which 
will be smeared over the segment. Discrete rings 
at segment ends in a region are. described at the 
end of all segment data for that region. 

E. Thermal Clue 41-44- A4 

This clue describes the type of’ thermal problem 

which exists in the segment. The- .user is. reminded . i. . 

that in a stability run only one of the two 
possible load cases can be thermal, and this clue 
pertains to that case. If there is ; no thermal 
load on the segment., the clue to be used is NOTH . 

If the thermal loading on the segment is of general, 
standard type, that is if there is variation 
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Column 


Format 


of temperature through the thickness, as. well as. 
in the coordinate directions, the clue- to.be •. ■ ., - 

used is THST . If the thermal load is such. that . , 

the variation is all in the coordinate, directions, 
and there is no thermal variation through the / 
thickness, the clue, to be used: is THCN . The last 
clue concerns a shell which is inhomogeneous in 
the meridional direction. This' is not really a 
thermal problem at all, but merely a manipulation . 
of the material property tables.: -If a structure 
has a wide variation in material properties in the : 
meridional direction, without this last option one 
must take short segments of constant properties : ' 
for analysis. With this option, however, the 
property variation is placed in the material 
property table, and expressed on the segment as v - 
a function of temperature. No thermal loads are - 
calculated, however, and the temperatures are--'' 1 
only used to interpolate for -material properties 
as integration is progressing' along the. segment-.' ■ 1 
Thus continual variation of properties in the 
meridional direction- is accommodated. - The clue ' 
for this option is THIN. : - 


F. Stress-free Temperature 5 

The value of the temperature (usually room tern- -..i • 

perature) at which the segment has no thermal 
stresses or distortions induced, is provided 
here. This is the temperature, at which the 
shell was manufactured. If. there is to. be.no : c 
thermal analysis, this value is not used and can r r v 
be set to zero (0.0). 




E10.1 
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Column 


G. Table control - Number of points in each of the 71-72 
following tables. 

This can vary from 2 to 30 depending upon the 
shell geometry' and loading. For a linearly : 
varying geometry and/or loading only 2 input 
points would be required. These- two' points 
would be the end points. For more general 
loading and/or geometry a larger number of 
points are required. Ih. particular, each 
abrupt change is specified by two points. • 

One should use as many points as necessary 
(up to 30) in order, to completely describe 
the problem, rather than using very short 
segments. 

5. Table of <J> or s Values 

A. Initial, intermediate and final values of 

<t> or s. Each point requires l4 . columns on a 
card and thus there can be 5 values per card 
and up to 6 cards to make, a total of up to 
30 points. The points to be specified are 
the beginning point of the segment, any point 
of discontinuity, and the end point of the 
segment. The input must be consistent with 
item G of the previous card. 

6. Table of Wall Crossection Geometry 

The contents of these cards (up to 6 cards per 
item below) are dependent upon the clues reg- :• 
istered on the Master Clue Card. If the shell 
to be described contains no reinforcing, the 
pertinent clue is item 4c, the Sheet Clue. The 
geometry is input and the stiffnesses are cal- 
culated by the program (see F-igure 2-10). The 
input is presented below as a function- of the 
Sheet Clue. 

- If the Sheet Clue is SING (single sheet-con-- --- r- 

struction) : 

A. Initial, intermediate and final values of wall 
thickness (h. ) at points defined by. table of ,<p 
or s values. 


Format 


12 


5E1^.7 


5E14.7 
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Figure 2-10 Shell Section Properties 












Column 


Format 


If the Sheet Clue is EQUA (equal-size face sheet 
shndwich) : 

A. Initial, intermediate and final values of face 
sheet thickness (h.) at points defined by table 
of f or s values. 

B. Initial, intermediate and final values of core 
thickness (t) at points defined by table of <j> 
or s values. 

If the Sheet Clue is UNEQ, (unequal-size face sheet 
sandwich) :• ' 

A. Initial, intermediate and final values of inner 
face sheet thickness (h. ) at points defined by 
table of <f> or s values. 

B. Initial, intermediate arid final values of core : 
thickness (t) at points defined by table of <f> 
or. s values. 

C. ' Initial, intermediate and final values of outer 

face sheet thickness (h ) at points defined by 
table of <J> or s values. 

If the shell is reinforced, the Sheet Clue will be 
BLAN. ' In this case it is the following, or Reinforce- 
ment Clue (item 4 d) which will determine the contents 
ofi card series 6. For- the reinforcement cases the 
geometry can be complex and varied, since all types 
of reinforcing are to be included. The reinforced 
shell-input is presented below as a function of the 
Reinforcement Clue. 

If the Reinforcement Clue is ST11 (single sheet 
reinforced by rings -and/or stringers) : - 

A. Initial, intermediate and final values of the 
torsional stiffness in the <}> direction (GJ ) at 
points defined by table of (}> or s values. ^ 

B. i Initial, intermediate and final values of the 

torsional stiffness in .the 9 direction (GJq) at 
points defined by table of f or s values. 

C. Initial, intermediate and final values of 

stringer area (A ) at points defined by table of 
<J> or s values. 9 " ' 

D. Initial, intermediate and final values of ring- 

area (A 0 ) at points defined by table of or s 
values. . • 

E. Initial, intermediate and final values of 

. stringer eccentricity (measured inwards from 
base shell centroid as positive) at points 
defined by table of 4> or s values. 
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5E14.7 
5El4 . 7 

5E14.7 

5E14.7 

5E1^=7 


5E14.7 

5Ell . 7 
5E14 . 7 
5E14 . 7 
5E14 . 7 



Column 

F. Initial, intermediate and final values of ring 

eccentricity (measured inwards from base shell 
centroid as positive) at points defined by ' 

table of or s values. - ! 

G. Initial, intermediate and final values of 
stringer moment of inertia (about base shell 
centroidal axis) at points defined by table 
of ♦ or s values. 

H. Initial, intermediate and final values of ring 
moment of inertia (about base shell centroidal 
axis) at points defined by table of (j> or s 
values. 

I. Initial, intermediate and final values of 
stringer spacing at points defined by table 

of tf) or s values. (Do not set to zero if no stringers.) 

J. Initial, intermediate and final values of. ring 
spacing at points defined by table of <f> or s 
values. (Do not set to zero if no ringa) 

K. Initial, intermediate and final values of base 
shell wall thickness (h.) at points defined by 
table of <p or s values. 

If the Reinforcement Clue is ST12 (equal face sheet 

sandwich reinforced by rings and/or stringers): 

A. through J. The items contained on these cards are 10 sets of 
those described for the ST11 clue above. 

K. Initial, intermediate and final values of base 
shell face sheet thickness (h. ) at points de- 
fined by table of <j> or s values. 

L. • Initial, intermediate and final values of base 

shell core thickness (t) at points defined by 
table of <j> or s values. 

If the Reinforcement Clue is ST13 (unequal face 

sheet sandwich reinforced by rings and/or stringers): 

A. through J. The items contained on these cards are 10 sets of 
those described for the ST11 clue above. 

K. Initial, intermediate and final values of base 
shell inner face sheet thickness (h. ) at points 
defined by table of <j> or s values. 

L. Initial, intermediate and final values of base 
shell core thickness (t) at points defined by 
table of ij> or s values. 
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Format 
5E14 . 7 

5E14 . 7 

5El4 . 7 

5El4 . 7 
5E14 . 7 
5El4 . 7 

5El4 . 7 
5El4 . 7 

5El4 . 7 

5El4 . 7 
5El4 . 7 

5El4 . 7 



Column 


Format 


M. Initial, intermediate and final values of base 
shell outer face sheet thickness (h ) at points 
defined by table of $ or s values. 

If the Reinforcement Clue is RWA1 (single sheet 
reinforced by a waffle rotated at an arbitrary angle 
from the meridional direction): 

A. Initial, intermediate and final values of waffle 
grid area at points defined by table of <j> or s 
values. 

B. Initial, intermediate and final values of waffle 
grid eccentricity (measured inwards from base 
shell centroid as positive) at points defined 
by table of <f> or s values. 

C. Initial, intermediate and final values of waffle 
grid moment of inertia (about base shell centroidal 
axis) at points defined by table of <j> or s values. 

D. Initial, intermediate and final values of waffle 
grid spacing at points defined by table of $ or 
s values. 

E. Initial, intermediate and final values of waffle 
grid rotation angle, p, (in radians from the 
meridional direction) at points defined by table 
of I or s values. 

F. Initial, intermediate and final values of extreme 
distance from base shell centroid to waffle outer 
edge (± value, positive inwards), at points 
defined by table of (p or s values. 

G • Initial, intermediate and final values of base 
shell wall thickness (h. ) at points defined by 
table of 4> or s values. 

If the Reinforcing Clue is RWA2 (equal face sheet 
' sandwich reinforced by a waffle rotated at an 
arbitrary angle from the meridional direction) : 

A. through F. The items contained on these cards are 6 sets of 
those described for the RWA1 clue above. 

G» and H. The items contained on these cards are 

those described for the ST12 clue above as items 
K. and L. 

If the Reinforcing Clue is RWA3 (unequal face sheet 
sandwich reinforced by a waffle rotated at an 
arbitrary angle from the meridional direction): 

A. through F. The items contained on these cards are 
those described for the RWA1 clue above. 

G. through I. The items contained on these cards are 3 sets of 

those described for the ST13 clue above as items 
K. through M. 


2 sets of 


6 sets of 


5El4 . 7 

5El4 . 7 
5El4 . 7 

5El4 . 7 
5El4 . 7 
5El4 . 7 

5E14.7 

5E14 . 7 

5El4 . 7 
5El4 . 7 

5El4 . 7 
5El4 . 7 
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Column 


Format 


If the Reinforcing Clue is ISG1 (single sheet 
reinforced by a general angle isogrid construction) : 

A. through D. The items contained on these cards 4 sets of 

are identical to those described for the RWA1 . 
clue above, but with reference to the isogrid. 

E. Initial, intermediate and final' values of the . 
isogrid angle, £, (in radians from the 
meridional direction, see Reference 1 
Appendix A). For the normal isogrid, the 

angle is 30° (inputin radians).. • ■ ■ r 

F. Initial, intermediate and final values of base- 

shell wall thickness (h. ) at points defined by : .. 
table of <|> or s values. ...... . 

If the Reinforcing Clue is.ISG2 (equal face sheet 
sandwich reinforced by a general angle isogrid 
construction)': 

A. through E. The items contained on these cards , 5 :?ets of 

are those described for the ISG1 clue above. 

F. and G. The items contained on .these cards are- 2- sets. of 

those described for the ST12 clue above as ' 

items K. and L. , ' . , 

If the Reinforcing Clue is ISG3 (unequal face sheet 
sandwich reinforced by a general angle isogrid 
construction): 

A. through E. The items contained on these cards • 5 sets of 

are those described for the ISG1 clue above. " ■ . - 


F. through H. The items contained. on these cards ' 3 sets of 

are those described for the ST13 clue above • • 
as items K. through M. . , 


If the Reinforcing Clue is STIC the following Hooke's. 
Laws will be used by the program for the description . 
of the shell wall: 

N e = K n e 0 + ^i 2 e b " c n k e n tr - 


K 22 e 4' + K 12 G 0 ~ C 22 k <j) _ N T4>- 

Y T o o ■ 


N = N = K y 
*0 0<t» ■ 33 


242 


5E14 . 7 
5E14 . 7 

5E14 . 7 

5E14 . 7 
5E14 . 7 

5El4 . 7 
5E14 . 7 



Column 


Format 


M e D n k e + D i2 k <i> + ( "n £ e ” 

o 

M <j) = D 22 k <t> + D 12 k 0 + C 22 £ <t> o ' **T4> 



-2D 


33 k 4>0 


Therefore the input is (see Ref. 1 Appendix A): 

A. Initial, intermediate and final values of ^ 
at points defined by table of <ji or s values. 

B. Initial, intermediate and final values of 

at points defined by table of if or s values. 

C. Initial, intermediate and final values of 

at points defined by table of <j> or s values! 

D. Initial, intermediate and final values of K-.-, 

•at points defined by table of <p or s -values. 

E. Initial, intermediate and final values of D^ 
at points defined by table of <p or s values! 
(Should be input as negative for sign convention. ) 

F. Initial, intermediate and final values of D^ 
at points defined by table of if or s values. 
(Should be input as negative for sign convention.) 

G. Initial, intermediate and final values of D 
at points defined by table of f or s values! 
(Should be input as negative for sign convention. ) 

H. Initial, intermediate and final values of D.,., 
at points defined by table of if or s values. 

I. Initial, intermediate and final values of 
at points defined by table of <p or s values. 

J. Initial, intermediate and final values of C ^ 
at points defined by table of <f> or s values! 

If the Reinforcing Clue is RWAF the following Hooke's 
Laws will be used by the program for the description 
of the shell wall: 


w e = K n e e + K i 2 e 4 > " c n k 0 " c i5 k <t> ” n ts 

- o o 


N <J> - *^22 £ <J> + K 12 £ 0 ' C 15 k 0 " C 22 %• _ N T<t> 

T T 


N 4>e ” K 33 r< P& 0 ~ 2C i 6 k <j>0 


5E1 1 *-.? 
5E14- . 7 
5E14- . 7 
5E1^.7 
5E1^.7 

5E14.7 

5E14.7 

5El4 . 7 
5E14 . 7 
5El4 . 7 
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Column 


Format 


M 0 .... \ ± k e + D 12 k <}> + C 11 e 8 + C 15 E <j> “ ^0 

o o 

M <J> = D 22 k <j> + D 12 k e + C 15- e 6 0 + C 22 G 4> o " • 

M <)>0 = “ 2D 33 N>e + C 1 6 y 4>0 q 

Therefore the input is (see Ref. 1 Appendix A).: 

A. through J. ' The items contained on these cards are 10 sets of 
those described for the ST10 clue above. 

K. Initial, intermediate and final values of C at 
points defined by table of <j> or s values. ^ 


L. Initial, intermediate and final values of C^g. at 
points defined by table of <)> or s values. 

7. Table of Input Stress Resultants ( Stability Only ) 

This card set is included only if the Prestress 
Input Clue (item I page 2 -15) is set to unity. 

A. Initial, intermediate and. final values of the 
hoop stress resultant, Nq , at points defined by 
table of cp or s values. 

B. Initial, intermediate and final values of the 
meridional stress resultant, Ify, at points 
defined by table of cp or s values. 

If the number of loading conditions (item G page 
2- l4) is set to two, then the above two card sets 
A and B refer to the first (fixed magnitude) load 
condition. They should be then followed by two 
more card sets, C and D, identifying the assumed 
Nq, and N ^ respectively, for the second (eigen- 
value) load condition. Do not use normalized 
values for buckling load or stress resultant 
assumptions since the program equations are not 
normalized, and thus such values will be too 
small for proper calculations. 

7. Table of Mass Densities . ( Vibrations Only ) 

Dimensions = weight/volume/acceleration due to 
gravity. 

The contents of these cards (up to 6 cards per item 
below) are dependent upon the clues registered on 
the Master Clue Card. If the shell to be described 
contains no reinforcing, the pertinent clue is item 
i|C, the sheet clue. The input is presented below 
as a function of the Sheet Clue. 


2 -kb 


5E14.7 
5Elt . 7 
5E14.7 

5E14.7 

5EII+.7 



Column Format 


If the Sheet Clue is SING (single sheet construction): 


A. Initial, intermediate and final values of the wall 
mass density at points defined by table of cp or s 
values. 

5E14.7 

If the Sheet Clue is EQUA or UNEQ, (sandwich 
construction): 


A. Initial, intermediate and final values of the 
core mass density at points defined by table of 
cp or s values. 

5E14.7 

B. Initial, intermediate and final values of the face 
sheet mass density at points defined by table of 
cp or s values. 

5E14.7 

If the shell is reinforced, the Sheet Clue will be 
BLAN. In this case it is the following, or 
Reinforcement Clue (item to) which will determine 
the contents of card series 7- The reinforced shell 
input is presented below as a function of the 
Reinforcement Clue. 


If the Reinforcement Clue is ST11 (single sheet 
reinforced by rings and/or stringers): 


A. Initial, intermediate and final "values of the 
base shell mass density at points defined by 
table of cp or s values. 

5E14.7 

B. Initial, intermediate and final values of the 
stringer mass density at points defined by 
table of cp or s values. 

5E14.7 

C. Initial, intermediate and final values of the 
ring mass density at points defined by table 
of cp or s values. 

5E14.7 

If the Reinforcement Clue is ST12 or ST13 (sandwich 
reinforced by rings and/or stringers): 


A. Initial, intermediate and final values of the 

base shell core mass density at points defined by 
table of cp or s values. 

5E14.7 

B. Initial, intermediate and final values of the 
- base' shell face* sheet mass' density 'at points ' 
defined by table of cp or s values. 

5E14.7 

C. and D. The items contained on these cards are 
those described for the ST11 clue as items 
B. and C. 

2 sets of 5E14.7 


2-^5 



Column 


Format 


If the Reinforcement Clue is RWA 1 or-ISGl (single ' • ' 

sheet reinforced by a rotated waffle or isogrid):" 

A. Initial, intermediate 1 and final values of the base , 5E14.7 

shell mass density at points defined by table of 

q) or s values. 

B. Initial, intermediate and final values . of the / , 5E14.7 

reinforcement mass density at points defined by,/;. 

table of 9 or s values. 

If the Reinforcement Clue is RWA 2 , RWA3, ISG2 or ISG3 
(sandwich reinforced by a rotated waffle or isogrid): 

A. Initial, intermediate and final values of the base 5E14-.7 

shell core mass density at points defined by table 

of tp or s values.. ' ■ 

B. Initial, intermediate and final values of the base 5E14-.7 

shell face sheet mass density at points defined" by " ’ 

table of 9 or s values. 

C. Initial, intermediate and- final values of the 5E14.7 

reinforcement mass density at points defined by 

table of cp or s’ values . 

If the Reinforcing Clue is ST10 or RWAF the program 
does not have any crossection description available 
to it for calculations . Thus the required input in 
card set 7 for these clues is (see Refs. 1 and 5): 

A. The meridional distributed load on the crossection 5E14.7 

due to rotation at points defined by table of 'cp or 
s values. For example for a layered section this 
can be expressed as: 

n ' - ■ "3 ■, 

f^ = / p.tiFr.t. cos 9 • 

cp i i i • 

i=l - 


B. The normal distributed load on the crossection due 
to rotation at ppints defined by table of cp or s 
values. For example, for a layered, section this.-, 
can be expressed as: : 
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5E14-.7 



Column 


Format 


C. The circumferential distributed moment load on ••••";> - *'• 5E14.T 

the crossectiori due to rotation at points .• . •- 

defined by table of cp or s values. For example 

for a layered section this can be expressed as:. 

n 

m e = *y] p i c ^ r i t i^i cos -{ P 

i=l 

Note: The sign convention for the above loads is • • 

as shown in Figures 2-11. 

D. The crossection parameter a^ at points defined by 5E14.7 

table of cp or s values. . 

E. The crossection parameter a, at points defined by 5®l^-7 

table of cp or s values. 

; F. The crossection .parameter a 2 at points defined, by 5 E l^-7 

table of cp or.s values-. v .,, . 

The crossection parameters a* can be defined as follows: 

a PC^^C. ‘ . j = °> 1 > 2 

"nhickness 

where Q is measured from the centroid ( inwards 
positive). 

8 . Loading Clue Card (These cards and card sequence 9 
are omitted under the following conditions: 

Stability - if the prestress input clue «( item I 
page 2- 15) is unity. 

Vibrations - if the Dynamic Analysis description 
clue (item A page 2- 18) is FREV, CRSP or CRSR. ) 

The contents of this card are numerical clues which 
alert the program to the types of loads that exist 
on the segment. In addition, if the clue indicates 
that some load does not exist, the appropriate 
cards in series 9 which would ordinarily contain • 
the numerical values of this load are omitted from 
the sequence. .. -..i. 



The series of cards 8 and 9 are repeated for the 
number of loading conditions indicated on the 
Program Control Card (item G) up to a maximum of 
2 (except in Vibrations - see below). If for 
one of these loading conditions/ no load exists 
on the segment, then a blank-'- Loading Clue Card 
is inserted in the sequence at this point and 
the load cards are entirely omitted (see Ref. 6). 

Only in a stability analysis can there possibly be 
two sets of cards 8 and 9- In vibrations analysis 
with Dynamic analysis description ©lues of VPRE, 
PCRS or PCSR only the prestress load state is 
described here, the dynamic loads being calculated 
internally by the program. 


The appropriate clues are as follows : 

A. Thermal Clue 

If there are no thermal loads (Item 4 e Is 
NOTH) the clue number is zero (0). 

If there is a standard thermal variation through 
the thickness (Item UE is THST) the clue number 
is four (U). 

If the temperature is constant through the 
thickness (Item UE is THCN) or if the in- 
homogeneous option is used (item 4 e is THIN) 
the clue number is one (l). 

In a stability analysis only one loading 
condition can contain temperature. Thus in 
the case of two load conditions on the 
shell the above item must be zero on one of 
the clue cards 

B. Enter a zero (o) in column 2. 


C. Meridional Load Clue (f ) 

4 1 

If there are no meridional loads, then the 
clue number is zero (0). 


If there are meridional loads, then the clue 
number is one (l) . 


D. Normal Load Clue (f ) - 

If there are.no normal loads, then -the clue 
number is zero (0). 


If there are normal loads,. then the clue 
number is one (l). 



Column 

E. Circumferential Moment Load Clue (m ) 5 

0 

If there are no circumferential moment loads , 
then the clue number is zero ( 0 ). 

If there are circumferential moment loads , then 
the clue number is one (l). 

F. Enter a zero (0) in column 6 6 

G. Any alphameric information (load description) 7-70 

9. Table of Applied Loads (see Figures 2 -lla, b for 

sign convention). 

The appropriate card sequence is given below as a 
function of the Loading Clues on card 8 .. If the Thermal 
Clue is one (l): 

A. Initial, intermediate and final values of the 
temperature of the shell at points defined by 
table of <j> or s values. (These values will be 
used either for a thermal problem where there 
is no thermal variation through the thickness 
{Clue = THCN), or to calculate varying material 
properties along the shell for an inhomogeneous 
problem {Clue = THIN}.) 

If the Thermal Clue is four ( 4 ) : 

A. Initial, intermediate and final values of the 
temperature T^ at points defined by table of 
♦ or s values. (The subscripts "nm" indicate 
temperature location - see below. ) 


T 

ii Inside 


$5$ 

Face Sheet ] 



1C 

oc 

T 

H 

II 1 III 
oneyco 


[L 

t - 


m f 

ace Sheet / 


T oo ' ^tside ”7 


Format 

II 

II 

16A4 


5E14.7 


5 E 1 I +.7 
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Format 


B. Initial, intermediate and final values of the 
temperature T at points defined by table of 
<(> or s values? 

C. Initial, intermediate and final values of the 
temperature T qc at points defined by table of 
♦ or s values? 

D. Initial, intermediate and final values of the 
temperature T qq at points defined by table of 
<j> or s values? 

If the Thermal Clue is zero (0), the above cards are 

omitted.,. . ■' j 

: * j 

If the Meridional Load Clue is one (l) : •' 

E. Initial, intermediate and final values of the 
meridional loads f at points defined by table 
of <p or s values . ^ 

If the Meridional Load Clue is zero (0), cards E are 

omitted. 

If the Normal Load Clue is one (l): 

F. Initial, intermediate and final values of the 
normal loads f at points defined by table of 
4> or s values . 

If the Normal Load Clue is zero (0), cards F are omitted. 

If the Circumferential Moment Load Clue is one (l): 

G. Initial, intermediate and final values of the 
circumferential moment loads at points 
defined by table of $ or s values . 

If the Circumferential Moment Load Clue is zero (0), 

cards G are omitted. 


5E1U.7 

5E1U.7 


5E1U.7 


5E1U.7 


5E1U.7 


5E1U.7 
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A > A ,w = A , and v 

/, n C 

in meridional plane 




A lie 
4 > 


Figure 2 -11c. Shell Element Geometry and Displacements 
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I 

Column 

10. Segment Topology Cards 

A. Segment number 1_5 

Number of the segment under consideration. 

B. Joint (i) 6-10- 

Joint associated with i^* 1 end of the segment 
(TIC). 

C. Joint (j) ■ 11-15 

Joint associated with the end of the 
segment (STOP). 

Since within a region the segments are all 
singly connected, the segment joint numbers 
should be in adjacent numerical pairs. That 
is, if joint ( j ) is 6, joint (i) could only 
be 5 or 7* This is true only within a region. 

In addition, the initial joint of each 
region must be 1 in segment topology 
numbering, and, the final joint of each 
region must be the last (highest) number 
in the segment topology numbering (see 
Figure 2-12) . The coordinate <p or s 
increases from TIC to STOP, i to j . 

The user is again advised to see 
Figures 2-3 to 2-9* 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus in 1-80 

(See General Notes - Data Debugging) ■ 

INTRA-REGION DISCRETE RING CARDS 

These cards, if any exist (region introductory card, 
item 1C), are placed at the end of all the segment 
data for the region. They contain the following 
information for each ring (in groups of 4 cards): 

1. Ring Location and Property Card 

A. Segment joint number to which the ring is 1-2' 12 

attached. If there is also a radial 
discontinuity at this location necessitating 
a kinematic link, the joint number of the 
ring is the independent joint number of the 
link (see Figure 2-13a) . 


Format 

15 

15 

15 
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Format 


B. Ring extensional stiffness (EA) 3-l6 

C. Ring bending stiffness about centroidal y 17-30 

axis (El ). See Figure 2-l4. 

y 

D. Ring cross-bending stiffness about centroidal 31-44 

axes (El ) . . ■ 

xy 

E. Ring torsional stiffness (GJ) 45-58 

F. Ring bending stiffness - about centroidal 59-72 


x axis (El ) . 

x 

Ring Geometry Card (see Figure 2-l4) 


A. Ring thermal coefficient of .expansion . (cO . . .1-12 . ; 

B. Ring centroidal radius (r ). 13-24 

C. Distance (±) between ring centroid and shear 25-36 

center in the x direction (x ). 

c 

D. Distance (±) between ring centroid and shear 37-48 

center in the y direction (y ). 

E. Offset (±) of ring shear center from attached _ 49-60 

shell joint (see item 1A) in the x direction (x). 

F. Offset (±) of ring shear center from attached 61-72 


shell joint (see item 1A) in the y direction (y) . 

Note; The shear center (S ) is the origin of the 
x, y axes shown in Figure 2-14. 

Ring Thermal Description Card (See Figure 2-l4) 

A. Distance (±) between ring shear center and l-l4 

the closest (innermost) radially measured 

ring point x^ (in x = radial direction). 

B. Distance (±) between ring shear center and. the. . - -.15-28 
farthest (outermost) radially measured ring 

point, Xq (in x = radial direction). 

Note; Items A & B cannot both be input, as zero. 


' El4". 7 
E14.7 

E14.7 

. E14.7 
E14.7 

E12.5 

E12.5 

E12 . 5 
E12 . 5 
E12 . 5 
E12.5 

Ei4.7 

E14.7 
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z 



Links Necessary- 



Independent joint 
of kinematic link 
and therefore Ring joint 
number . 

Kinematic Link Necessary 


(a) Segment or Simple Region Rings 


kinematic link 1 



Independent joint for all 
.kinematic links and 
therefore Ring joint 
number. 


(b.) . Possible Region Multi- connected Rings 


Figure 2 - 13 . Discrete Ring Topology 






Column 


Format 


c. 

Innermost point temeperature (T^) 

29-k2 

Ell. 7 

D. 

Outermost point temperature- (T q) 

^3-56 

Ell. 7 

E. . 

Ring stress-free temperature (room or manufac- 
turing temperature (see Note for Master Clue 
Card item F). 

57-70 

Ell. 7 

F. 

Load Case Number (Stability analysis only) 

71-72 

12 


Enter the load case number (either 1 or 2) 
which contains the above described thermal 
load. This should be consistent with the 
segment thermal load order. 

Note : Thermal loads on the rings are input on 

the above cards. Mechanical loads on the ring 
must be input as a load distribution on the 
segments connecting to the ring. A higher"" 
intensity loading may be distributed near the 
segment ends connected to the ring. 

The number of sets of discrete ring cards 
must be equal the number specified in item 1C, 
of the region introductory card. 


Ring Prestress Card (Stability analysis only ) 


A. 

Radius of revolution at the segment joint 
which is used to identify the ring (item 1A). 

l-ll 

Ell. 7 

B. 

Hoop prestress in the ring due to load 
condition 1 . (Enter only if the prestress input 
clue, item I page 2-15, is set to unity.) 

15-28 

Ell. 7 

C. 

Hoop prestress in the ring due to load 

29-12 

Ell. 7 


condition 2. (Enter only if the prestress input 
clue, item I page 2 - 15 , is set to unity , and 
there are two loading conditions. ) 


Ring Mass Card (Vibrations analysis only ) 


A. 

Radius of revolution at the segment joint 
which is used to identify the ring (item 1A). 

l-ll 

Ell. 7 

B'. 

• Ring mass density (p) - - - 

- - 15-28 

- -Ell . 7 

C. 

Ring Young's Modulus (E) 

29-12 

Ell. 7 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D 
(See General Notes - Data Debugging) 

minus 

in I -80 
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Format 


INTRA-REGION KINEMATIC LINK CARDS 

These cards, if any exist (region introductory card, 

item IB), are placed at the end of all the discrete 

ring data for the region. They contain the following o> 

information: 

A. Joint (j) dependent joint • 1-2 12 

B. Joint (i) independent joint - > 3-^ 12 

For intra-region kinematic links these joints 
must he in consecutive decending order. That 
is, joint (j) should always he greater than 
joint (i) hy one. ' 

C. Angle y in radians (see : Figure 2-15) . - 5-19 El4- . 7 

y cannot equal 0 or it. 

The angle y describes the orientation of the 
link; it is the inclination angle of the link. -,■< 
from the vertical (Z axis). 

The number of kinematic link cards must 
equal the number specified in item IB, . 
of the region introductory card. 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus in 1-80 

(See General Notes - Data Debugging) 

REGION JOINT CONTROL DATA ' ' 

These cards are placed at the end of all the data for 
all regions ... 

1. Joint Control Data Card' . . 

A. Number of region joints 1-5 15 

Total number of region joints (Max. = 30). ’ 

B. Number of discrete rings 6-10 15 

Total number of discrete rings between . . 
regions. 

C. Number of kinematic links 11-15 15 

Total number of kinematic links between - 
regions. [ 
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Discrete Ring Cards (inter-region, if any exist) 
A. Ring Location and Property Card 

Column "■ ■ 

Format 

a. 

Region joint number to which the ring is 
attached. If there are also radial 
discontinuities at this location neces- 
sitating kinematic links, or multi - 
connections, the joint number of the ring 
is the independent joint number of ail 
the links (see Figure 2-13b). 

1-2 

12 

b. 

Ring extensional stiffness (EA) 

3-16 

E14.7 

c. 

Ring bending stiffness about centroidal 
y axis (EI^). See Figure 2-l4. 

17-30 

E14.7 

d. 

Ring cross-bending stiffness about 

centroidal axes (El ). 

v xy 

31-44 

E14.7 

e. 

Ring torsional stiffness (GJ) 

45-58 

E14.7 

f. Ring bending stiffness about centroidal x 
axis (EI X )* 

B. Ring Geometry Card (see Figure 2-l4) 

59-72 . . 

E14.7 

a. 

Ring thermal coefficient of expansion 
<«„). 

1-12 

E12.5 

b. 

Ring centroidal radius (r ) . 

• c 

13-24 

E12.5 

c. 

Distance (±) between ring centroid and 

shear center in the x direction (x ). 

c 

25-36 

E12.5 

d. 

Distance (±) between ring centroid and 
shear 1 center in the y direction (y c ). 

37-48 

' E12.5 

e. 

Offset (±) of ring shear center from 
attached shell joint (see item Aa) in the 
x direction ("x). 

49-60 

E12.5 

f. 

Offset (±) of ring shear center from 
shell joint (see item Aa) in the y 
direction (y). 

Note: The shear center (s) is the origin 

of the x,y axes shown in Figure 2-14. 

61-72 

‘ E12 ; 5 


2-61 



Column 


Format 


C. Ring Thermal Description Card (See Figure 2-14) j v 

a. Distance (±) between ring shear center l-l4 El4.7 

and the closest (innermost) radially 

measured ring point, Xj (in x = radial 
direction) . 

b. Distance (±) between ring shear center and 15^28 El4.7 

the farthest (outermost) radially measured 

ring point, xq (in x = radial direction). , , 

Note: Items a & b cannot both be input as zero. 

c. Innermost point temperature ■ (Tj )■. . 29-42 E14.7 

d. Outermost point temperature (Tq). 43-56 E14.7 

e. Ring stress-free temperature (room or 57-70 El4.7 

manufacuturing temperature (see Note for 

Master Clue Card item F). • 

f. Load Case Number (Stability analysis only ) 71-72 - 12 


Enter the load case number (either 1 or 2) 

* i . 

which contains the above described thermal 
load. This should be consistent with the 
segment thermal load order. 

Note : Thermal loads on the ring are input 
on the above cards ; Mechanical loads on . 

the ring must be input as line loads or ■ 
moments on the shell "joint (itemAa) 
assigned to the ring. 

The number of sets of discrete ring cards 
must equal the number specified in item IB, 
of the Joint Control Data Card. 

D. Ring Prestress Card (Stability analysis only ) 

a. Radius of revolution . at the region joint l-l4 El4.7 

which is used to identify the ring ; (ItemAa)'; ' 

b. Hoop prestress in the ring due to lpad 15?-28 El4.7 

condition 1. (Enter only if the prestress 

input clue, item I page 2-15, is se-fc to 

unity . ) . , t , . ’ 

c. Hoop prestress in the ring due to load ! ’ 29-42 1 1 El4.7 

condition 2. (Enter only if the prestress"- <• 

input clue, item I page 2 - 15 ,' is se,t ' to :S 

unity , and there are two loading conditions. ) 
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Column 


Format 



Rine Mass Card (Vibrations analysis onlv) 


a. 

Radius of revolution at the region 
joint which is used to identify the 

1-.14 


ring (item Aa). 


b. 

Ring mass density (p) 

15-28 

c. 

Ring Young's Modulus (E) ' • 

29-4-2 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D 
(See General Notes - Data Debugging) 

minus in 


Kinematic Link Cards (inter-region, if any exist) 


A. Joint ( j) dependent joint 

B. Joint (i) independent joint 


1-2 

3-4 


For kinematic links between regions there 
are no restrictions upon joint numbering. 

The only restriction is that between successive- . 
kinematic link data cards the. (j) joint entry 
should be in increasing order (not necessarily 
consecutive). For example: 



i 


Must be 
in 

increasing 

order. 




C. Angle 7 in radians- (see Figure 2-15) 


5-19 


7 cannot -equal 0 or it. 

The ahgle 7 describes "the orientation of the 
link; it is the inclination angle of the 
link .from the vertical (Z axis). The number 
of kinematic link cards must equal the number 
specified in item 1C of the s joint Control , 
Data Card. j 


ElU.7 

ei4.7 

ei4.7 

1-80 


12 

12 


E14.7 
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Column 


Format 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus in 1-80 

(See General Notes - Data Debugging) 

Boundary Condition Cards (Joint Data - one card 
per region joint) 

This is the first set of boundary conditions. 

If the number of boundary condition sets is 
input as greater than 1 (item H page 2- 14) , 
they will be used for the static analysis 
and for the n = 0 eigenvalue analysis. If 
the eigenvalue harmonic search is specified 
to start with a harmonic other than zero 
(item D page 2- 14) then they will be used . 
for static analysis only. If the number 
of boundary conditipns is input as unity, 
these will be used throughout the analysis. 


A. Joint Number 

B. Joint component conditions on: 

1-2 

12 

!) Ap 

3-4 

F2.0 

2) or (see Figure 2-16) 

5-6 

F2.0 

3b ^ or \ 

7-8 

F2.0 

» % : 

9-10 

F2.0 

C. Angle a in radians 

11-24 

E14.1 


To be used only in conjunction with a 2 or 
3 code. 


Note: There must be as many boundary con- 

dition cards as there are joints as indicated 
in item 1A of the Joint Control Data Card. 




i 

i 


0 


0 


0 


+1 


0 


0 


0 

0 


+1 


0 


! o 


/ 

•i— J 

E-< 

<1 

s 

| ■ <r r r i ) 


A z. 

1 


< 


(r .-r. ) cot y 
1 J i 


■ a r . 
1 

1 +1 


fi 0. 


Figure 2-15. Kinematic Link 
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Figure 2-l6b. ; Provision for Local Rotations 
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There are 4 different codes that are used 
to prescribe joint component conditions. 
They are: 

a. 0 = no displacement allowed. 

b. 1 = displacement allowed in the 

indicated direction. 

c. 2 = A^ and A^ are rotated through an 

angle or(ir/2 - cc). and become A^. 
and A. respectively, while a 
displacement is allowed in the 
direction. 


d. 3 = A_, and A^are rotated through an angle 
of (ir/2 - a) and become A^ and A 
respectively, while a displacement 
is allowed in the A^ direction. 

See Figure 2-l6 for a geometric • 
explanation of codes 2 and 3. 


When using rotation codes: 


Code 2 can exist only as A^ coding. 
Code 3 can exist only as A,, coding. 


Codes 0 and 1 can appear in either 
column 4 or column 6, in addition 
to columns 8 and 10. Thus, there 
are twelve possible boundary con- 
ditions when rotation codes are 
used , (a = <j> for table below.) 



Free edge (possible 
to apply shear and/ 
or membrane loads) 

Aq = 0, normal 
support (possible 
to apply membrane 
load)- - - 

Ajj = 0, membrane 
support (possible to 
apply shear load) 

1, 2,3,1 

1, 2,0,1 

1,0, 3,1 

0, 2,3,0 

0,2, 0,0 . 

0,0, 3,0 

0, 2,3,1 

0,2,0, 1 

0, 0,3,1 
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Column 


Format 


Apex boundary conditions : 

Since the closed apex ( 4> = 0 ) angle 
is not acceptable input , the apex 
boundary conditions must be simulated 
at a small tf> angle. These boundary 
conditions vary per Fourier harmonic,, 
and are as follows : 


n 


Az, 


. n e 

angle 

0.0 

0 

0 

3 

0 

a = <t> 

1.0 

1 

2 

0 

l 

Q 

ii 

-e- 

>2.0 

0 

0 

0 

0 

— 


For a further discussion see 
Reference 1 , Appendix C. 


General Notes : 

1) To establish a datum for measuring 
displacement, free body motion must be 
eliminated from the structure. This should 
be accomplished by suitably applied boundary 
conditions. 

2) The ability of a dependent joint in a 
kinematic link to prescribe motion independ- 
ently should be removed by setting all 
boundary conditions of that joint to zero. 
See pages 2-63- and 2-74. 

3) In the axisymmetric static and 
eigenvalue analysis all torsion degrees 
of freedom should be removed. A second 
boundary condition set, allowing torsional 
freedom, should be used for the other 
harmonic eigenvalue analyses. 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D 

• (See General Notes - Data Debugging) minus in 1-80 
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Column 

JOINT LOAD DATA (These card sets are omitted if the 
Prestress input clue is set to unity in a stability 
analysis, or if the Dynamic analysis description 
clue is FREV, CRSP or CRSR in a vibrations analysis. ) 

1. Load Control Data Card 

A. Number of Joint Loads 1-4- 

Tot al number of joint loads in analysis. 

(Line loads can only be applied to region 
joints.) 

B. Any alphameric information (load description) 5-68 

Note ; If there are no Joint Loads for the 
structure, card 1 of the JOINT LOAD DATA 
is blank and card set 2 is omitted. 

2. Joint Load Cards (as many as in item 1A above) 

A. Load Condition Number 1-5 

Number of loading condition in which the line 
• load exists. (See Program Control Carl item 

2G. For a vibrations analysis this number 
is always unity. ) 

B. Row Identification 6-10 

The identification is the location of the 
degree of freedom at which the load is 
applied. This is obtained by counting the 
non-zero codes entered in the Boundary 
Conditions Cards, starting with Joint 1; 

T, Z, R, Joint 2; T, Z, etc., 

and stopping at the joint and degree of 
freedom where the line load is to be 
. applied. The location number of this 
degree of freedom is the information 
necessary. 

C. Applied Joint Load 11-24- 

The input is 2irr times the running load 
in lb. /in. In tfie particular case of the 
axial axisymmetric load, this is simply 
the net force. For sign convention see 
Figures 2-11, 2-16 and 2-17- 

D-A-S-H S-E-P’-A-R-A-T-O-R C-A-R-D minus 

(See General Notes - Data Debugging) 


Format 


14 


16A4 


15 


15 


E14.7 


in 1-80 
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Format 


SUPPLEMENTARY BOUNDARY CONDITION, DATA (The maximum 
number of sets of this data can be two. If the 
number of boundary conditions is input -as 2 (item H 
page 2- 14 ) the set of boundary conditions here will 
be used on all the eigenvalue analyses other than 
n=0. If the number of boundary conditions is input as 
3 (item H page 2- i4 ) the first set of the boundary 
conditions here will be used on the first n f 0 
eigenvalue analysis, while the second set here will be 
used on all the following, eigenvalue analyses.. 

1. Boundary Condition Card Set 2 (OMtt unless the 
number of boundary conditions is iput as 2 2) 

Joint Data - one card per region joint (See 
notes on page 2- 64). 


A. 

Joint Number 

• 1-2 

12 

B. 

Joint component conditions on: 




1 ) a t 

3-4 

F2.0 


2) Az or A^ .( see Figure 2-16) 

5-6 

F2.0 


3) A r or A q 

7-8 

F2.0 


4) be 

9-10 

F2.0 

C. 

Angle ot in radians 

11-24 

E14.1 


Note : There must be as many .boundary . condition ' • • •’ 

cards as there are joints as indicated in item >. 

1A of the Joint Control Data Card. 

... . . > • 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D . minus in 1-80 

(See General Notes - Data Debugging) 

2. Boundary Condition Card Set 3 (Omitt unless the ; , 

number of boundary conditions is input as 3) 

Joint Data - one card per region joint (see notes on 
page 2- 64). 

A. Joint Number . 1-2 12 

B. Joint component condition on: 

1) A t 

2) Az or Ajj(see Figure 2-16) 

3) A r or A q 

4) fig. . ------- - ... ; 

C. Angle a in radians : 

Note : There must be as many boundaiy condition 
cards as there are joints as indicated in item 
1A of the Joint Control Data Card. 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D -• minus in 1-80 

(See General Notes - Data Debugging) 
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3-4 F2.0 

.. .5.-6 F2.0 

7-8 F2.0 

------ 9-10 - - ■ - - J 2-.0 - 

11-24 E14.1 



Axial Load (Z) 



Program Input Load = 2nr x (F ) 

— — O L 


lb/in 


(negative in 
pictured direction) 


Z 




(f t ) 


Program Input Load = 2iTr o x (M) 
»M_ = M in-lb/in 


Figure 2-17a. Line Loading for Harmonic n = 0.0 
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Figure 2-17b Line Loading for Harmonic n = 1.0 
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It is hoped that the user is now able to use the STARS-2B,-2V programs td good 
advantage. They are : powerful toolsi, which will increase in value to the user 
as he uses them. One of the more complex . are as of usage is the description 
of topology, especially when involved with rotation codes and joint loads. 

An illustrative example of a Y joint representation is therefore presented 
below. (See Figures 2-18 and 2-19 for the structure and idealization.) 

The idealized structure contains four regions and two kinematic links. The 
joints are numbered from 1 to 7* Membrane loading is applied to joint's 1 
and 5 and the structure is supported by membrane action at joint 7- All. 
regions must be coupled* 

The second card in each region description (topology card) is as follows: 


Region 

Joint (i) 

Joint (j) 

1 

1 

2 

2 

2 

3 

3 

• '' 4 

5 

4; 

6 ' 4 

7 


The Joint Control Data card would contain a, 7 in column 5 and a 2 in column 15. 

In this example, the restraints at joints 1, ^5 and 7 must be rotated from the 
fixed (global system) to a local system such that membr.ane action may be 
applied. In addition, joints 3 and 4, and 6 and 3 are to be coupled with 

kinematic links. Thus, the motions of joints 4 and 6 are dependent upon the 

motion of joint 3- This dependence will be insured by using 2 kinematic link 

cards and setting the displacements of joints 4 and 6 equal to zero. It should 

be noted in this particular case that the motion of joints. 4 and 6 is not being 
equated to zero, but rather, the ability to prescribe motion independently is 
being removed. The required data has the following appearance. 

A. Kinematic Link Cards (2 cards) 


Joint (j) 

4 

6 


Joint (i) Angle 

3 ■ 7+tt 

3 7 


Note : In a double link of the 

type shown in Figure 2-19, one 
joint must be consistently 
independent (joint 3 in example). 






B. Boundary Condition Cards for higher harmonic eigenvalue analysis 
( 7 cards ) 

Joint T Z _R_ 9_ Angle 

1 1 2 3 1 O-l 

2 1.111 

3 1111 

k 0 0 0 0 

5 1 2 3 1 Og 

6 0 0 0 0 

7 0 0 3 1 a 3 

The Load Control Data card would contain a 2 in column k. 

The external membrane loads (one load problem assumed) are applied to the 
structure through the Joint Load Cards which, in this example, would appear 
as (2 cards): 

Row Load 

2 ^"o x Membrane Load (lb. /in.) 

1^ ^ 7Tr 0 x Membrane Load (lb./in.) 


Problem No. 
1 
1 
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SECTION 3 


OUTPUT INFORMATION 


The output of the STARS -2B, -2V programs is straightforward, however a 
description is in order since the user should learn the significance 
of the various checks that are provided. It is important to point 
out that the output of the program will include a print-out of the 
input data. This gives the user the opportunity to check whether 
or not the input data was correct. In the detailed description of 
the complete output which follows, the user should refer to the output 
of the problems in Section 4 as an example. 


The title page of the output contains all the data from the General 
Introductory Cards, prominently placed, and needs no comment. The 
next page of the output contains the first region Identification 
Card in the center.. The following output is then presented for each 
segment in this region (in order of appearance): 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 
9. 

10 . 

11 . 

12 . 

13. 

14. 


Contents of segment Identification Card. 

Contents of MAGIC Control Card 

Contents of Geometric Description Card (Cards) 

Contents of Master Clue Card 

The material property table used for the segment 
Crossection description table 
Temperature load table (if any) 

Distributed load tables (printed per loading condition; if any) 
Prestress input (if any) STABILITY ONIY 
Segment influence coefficients (MAGIC output) 

Segment stiffness matrix 
Stiffness matrix symmetry check 
Segment load matrices 

Radius of revolution at .1 and j ends of segment. 
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Item 12, the stiffness matrix symmetry check, is a check upon the 
validity of segment sizing and: the accumulation of round-off error. 

For perfect symmetry to exist, it is necessary to have zeros above 
the main diagonal, and zeros or ones below the diagonal. The amount 
of error induced by improper sizing or round-off is related to the 
amount that the off-diagonal terms in the lower triangle differ 
from unity. An attempt should be made to keep the upper limit on 
this difference at one percent (maximum number in lower triangle 
should be 0.1010... E Ol) . 

As mentioned previously, items 1-14 are repeated for all the segments 
within region one. The radius of revolution at every joint should be 
checked at corresponding joints of adjacent segments to make sure 
that proper coupling has been specified. At this stage in the output, 
the topology of the segments within the region, and the description 
of the intra-region rings and kinematic links is presented. 

Next the region matrices are presented. Given in order are: the 

region stiffness matrix, the stiffness matrix symmetry check, and 
the region load matrices. Again the numerical round-off, evident 
in the symmetry check, should be kept to a maximum of one percent 
(O.IOIO... E 01 in lower triangle). The output to this point, that 
is, sets of items l-l4, segment topology and links, and the region 
matrices, are now repeated as a group for each region within the 
structure. When this is completed, the region topology is presented. 

The next items to be provided by the output are the descriptions of 
the inter-region rings and kinematic links and the boundary conditions. 
At external points of the structure these are physical boundary 
conditions. At internal points they merely state the fact that no 
restraint exists and the joint in question is free to move. The last 
column in this set, gives the angle a, which is zero unless a rotation 
code is indicated . It is important to refer to Figure 2-6 once more and 
point out that a represents a rotation of the coordinate system. 
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There are a variety of errors that can be made in submitting input data 
The STARS -2B,-2V programs are set so that as an error is detected, (input 
or in ythe matrix .calculations) the. program is stopped. Therefore to avoid 
delays... in getting an answer, the SATELLITE - IB , - IV programs were created to 
debug .input data, and to check for other possible inconsistencies. The use 
of these programs is described in Reference 3, and all the STARS error messages 
are contained in the, listings therein. 

After having corrected any errors, so that it is nc>w possible for the program 
to run to completion, the problem solution can be discussed. From the above 
point in the output, further output will depend upon the intermediate print 
clue. If this clue is set for full output all stiffness and prestressed (or 
dynamic in the case of vibrations) stiffness matrices will be printed for all 
segments, regions, and structures in all harmonics and all estimates. If the 
print clue is set for abbreviated print, only that noted below will be printed: 

1. The applied line loads. 

2. The region end deflections per problem, in numerical order, starting 
from region joint 1. 

3. A description of the calculated prestress state (if any). This will 
include. the hoop and meridional stress resultants and the meridional 
rotation; The print will be for, two loading cases if so input. 

b. New boundary conditions for eigenvalue analyses (if any). 

5. The current harmonic number and the current estimate (iteration) number 
for that harmonic. 

6. All the. eigenvalues for the estimated prestress or frequency state, and 
: the current load (frequency) multiplier for the next estimate. 

7. Prints 5 and 6 are repeated until convergence (or a cutoff value) is 
reached, for . this harmonic. At this point the load multiplier represents 
the lowest eigenvalue, (The frequency is printed in a vibrations analysis.) 

- - . Two- eigenvectors- are. then printed- in global coordinates, and expanded- in - 

local deformation coordinates throughout all the segments in shell. In a 
stability analysis these two vectors will correspond to the lowest two 
... eigenvalues. In a vibrations analysis they will correspond to the two 
lowest eigenvalues, or the two eigenvalues closest to the estimated 
frequency, as requested by input. 
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8. Prints k~J are repeated for all the harmonics under investigation. 
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SECTION 4 


EXAMPLES OF INPUT AND SOLUTIONS 


This section contains the full input and output of two example test problems 
for the STARS -2B program, and one example test problem for the STARS-2V 
program. Since the test problem for the STARS-2S program C6) exemplified 
all crossection options available in the programs, the present test problems 
will be more physical in nature, while still exemplifying various analysis 
options. A comparative description of the present test problems with other 
analyses is presented in Reference 1. 

The first stability test problem is presented in Figure 4-1. A prolate 
spheroid is loaded in hydrostatic compression and the stability load is being 
sought. Since it is known Cl] that the critical harmonic for this case is 
n = 3 , only the 2nd and 3rd harmonics are investigated. The input is set 
so that the program will calculate the prestress state, and a cutoff value is 
set at 1.5. No intermediate output is required. It should be noted in this 
problem, that the boundary conditions to simulate an apex must be changed for 
harmonics n = 2, 3 from those used for the n = 0 static analysis (see Section 2) 
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ELLIPSE stability test 
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The second stability test problem is presented in Figure 4-2. Actually the 
test problem, involves only the panel buckling mode of the cylindrical shell, 
whereas the full stability analysis is discussed in Reference 1. The 
reinforced (smeared stringers) panel is loaded with a constant 31 psi. internal 
pressure, and a variable compressive end load to buckling. The input is set for 
the program to receive the prestress states from cards. There are two loading 
conditions: the constant pressure load , and the variable end compression. 

Since no axisymmetric static analysis is necessary, only one set of boundary 
conditions is sufficient. The input is set to include full intermediate output. 
Since it is known W that the critical harmonic is n=10, the search is set 
between n=9 and n=ll. 
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Notes: 

Loading 

1) Compressive end load: N 

2) Internal stabilizing pressure: 31 psi 
Boundary Conditions 

1) Ends are simply supported 


Figure 4-2 Test Problem 2 
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segment number I segment CODE Si SHUTTLE cylinder 



matrix x and y itr»nsposcd) hxgic- output 
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THE REDUCED STIFFNESS MATRIX 
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RES ION NUMBER 3 
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there aRe « segments and o kinehatjc links "ithin this Re6ion 
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the reduced stiffness matrix 
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there *Re i Segments and o kinematic links within this region 
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there *re i segments aNo o kinematic links »ithin this region 
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there ARE I SEGMENTS and 0 KINEMATIC links *|ThIN this region 








there are i segments ano a kinematic links "ithin this region 
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THE REDUCED stiffness MATRIX 
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THE REOUCEO STIFFNESS HaTRIX 
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there are > segments and o kinematic links *ithIn this region 
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there are i segments and o kinematic links «iThin this Region 
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REOUCED STIFFNESS MATRIX 
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The third test problem involves the investigation of vibration frequency 
characteristics of a prestressed spherical cap as shown in Figure 4-3. In 
this problem five harmonics are searched (n = 0-4) and thus three sets of 
boundary conditions are utilized. The input is set so that the two lowest 
eigenvectors are investigated in each harmonic. No intermediate print is 
required. The prestress load is a uniform external pressure of p = 1466 psi. 

It is interesting to note in this problem that the prestress pressure is 
rather close to the buckling load for the n=2 harmonic. The buckling load of a 
shell can be found using a vibrations analysis, by finding the loading which will 
cause a zero frequency. Further analytical discussions and comparisons for this 
problem are available in Reference 1. 


4-258 




Figure 4-3 Test Problem 3 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS. TO SI UNITS 

' , - , / 

The International System of Units (SI) was adopted by the 

Eleventh General Conference on Weights and Measures in I960. Conversion 
! factors for the units used in this report are given in the following 
table: 


Physical quantity 

U.S. Customary 
Unit 

mPmSM&wM 

SI Unit 
(**) 

Length 

in. - 

0^0254 

meters (m) 

Stress modulus 

ksi 

6.895 x 10 6 

new tons /mete 
(N/m^) 

Stress resultant 

lbf/ in. 

! 

175-1 

newtons/mete 

(N/m) 

Temperature change 

°F 

5/9 

Kelvin (K) 


* Multiply value given in U.S. Customary Unit by conversion factor to 
obtain , equivalent value in SI Units. 

** Prefixes to indicate multiple of units are as follows: \ 


Prefix 

Multiple 

giga (G) 

10 9 

mega (M) 

10 6 

kilo (k) 

10 3 

deci (d) 

10" 1 

centi (c) 

10- 2 

milli (m) 
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